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ABSTRACT 



The low-speed bluff body studies conducted on aeroelastic 
models were extended by redesign of the installation in the West 
Coast Research Corporation^ subsonic wind tunnel at the U. S. 
Naval Postgraduate School, Monterey, California* Improvements 
included stiffer support base, redesign of the spring constraint, 
and the use of eddy-current type displacement gages. In addition, 
the model* s lateral response data were monitored, conditioned and 
recorded on one-inch magnetic tape in a voltage analog form. 

The analog data were converted to digital form at a sampling 
rate of 500 numbers per second for a period of eight seconds per 
case in order to permit subsequent digital analysis in the Control 
Data Corporation 1604 computer. The results of this digital 
computation provided a numerical estimate of the spectral 
distribution, which gives an insight into the character of the 
random response of the aeroelastic structure due to the 
aerodynamically induced unsteady air load. 
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1 . 



Introduction . 



The use of aeroelastic type models has been widespread and is 
well known for estimating flutter boundaries on complete aircraft. 
Early airplane flutter models were primarily concerned with matching 
vibration frequencies, mode shapes, and mass scaling without 
concern for compressibility or viscosity effects. The advent of 
supersonic aircraft brought on an impetus to scale Mach number, 
hence the design and testing techniques on aeroelastic models 
changed in accord with the usage of transonic and supersonic 
wind tunnels [ 1 ] . The fairly recent importance of ground wind 
loads on missiles has brought emphasis upon bluff-body type 
aeroelastic models with similar scaling factors, except that 
Reynolds number is substituted for Mach number. Hence, a revision 
occurs in testing missile elastic models since pressurized low- 
speed wind tunnels are important in order to provide proper 
Reynolds number matching. 

In addition to studying specific three-dimensional bluff bodies, 
a need exists for determining and evaluating both the forcing 
functions and structural response characteristics of idealized 
single degree of freedom bluff body aeroelastic models with cross- 
sections as a test parameter. Consistent with this type of work, 
a program started in 1965 by Breckon and Hite [10] has been extended 
and the test installation and technique has been upgraded in order 
to provide a better degree of insight into bluff body model response. 
Engineering -type improvements included a stiffer support base, 
redesign of the spring constraint, and the use of eddy-current type 
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displacement gages. These gages operate as variable inductance 

instruments where inductance changes are a function of air gap. 

& 

Specifically, with regard to the antenna project the horizontal 
members, termed boom boards, have exhibited undesirable wind induced 
oscillations. Ideally the boom board structure could be visualized 
as a constant cross-section beam having some unknown end fixity 
(approximately a simple support) with an evenly spaced elastic load 
over its length and acted upon by an unknown random aerodynamic 
forcing function resulting from steady wind or gusts. 

The boom boards have cross sections of circular cylinders, 
squares, and rectangles. They all have the common feature of being 
bluff bodies, which means that the flow does not progress smoothly 
over the body, but separates from the surface in an unsteady manner 
which gives rise to the generation of fluctuating air loads. RMS 
amplitude of the response increases with increased wind velocity. 

The lightly damped structure acts as a band-pass filter to the 
random forcing function [2] . Usually the structure responds only 
to the forcing function centered about the resonant frequency of 
the structure. The principle response is essentially at the body’s 
natural frequency . 

The response of the model differs from that of the usual 
flutter model since bluff body response is associated with random 
type signals, while airplane flutter models involve harmonic type 
outputs. The analysis of random signals can be approached by pure 
analog signal processing methods where the random signal is 
recorded on magnetic tape and then played back through a harmonic 
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analyzer. Drawbacks to this technique have usually been that high 
quality filter circuits are needed when operating in the two cycle 
per second band width range; the lower limit on frequencies is 
between 20 and 30 cycles per second, and only power spectral density 
results are obtained. The auto-correlation function has not yet 
been obtained very satisfactorily by pure analog means. 

Another technique used in handling random data is the hybrid 
approach; this method was utilized in this survey. The data were 
recorded in analog voltage form and then converted to digital form 
utilizing a sampling rate consistent with the frequency content of 
the signals. A sample rate of 300 numbers per second allows 
spectral distributions up to 230 cycles per second, while a 
digital record of eight seconds in duration provides a lower bound 
on spectra of about one cycle per second. The conversion to 
digital form utilized an analog/digital converter in conjunction 
with a Control Data Corporation (C.D.C.) 160 computer. Once in 
digital form, the binary tapes were processed through a C f D.C. 1604 
computer. The output results were estimates of auto-correlation 
and power spectral densities. 

The low-speed wind tunnel used for these investigations was 
readjusted and calibrated during the summer of 1963 and the flow 
traits were improved from those reported by Breckon and Hite [ 10] . 
Main changes included removing the flow swirl and inclination in 
the test section by means of corner turning vane adjustments. A 
five orifice yaw head was designed and constructed for this purpose 
and allowed measurements of flow angularity to +0.02 degrees. 
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Turbulence measurements were made in February, 1966 using a 
0.0002 inch diameter tungsten hot wire probe. In addition to the 
usual RMS measurements, the turbulence measurements were processed 
through the analog to digital data handling system in order to 
provide a measure of the scale of turbulence and the spectral 
distribution of random flow eddies. Although the addition of fine 
mesh screens in the tunnel settling section are desirable, the 
tunnel in its present condition is quite usable for measurements 
of a character like those reported herein. 

The author worked very closely with Lieutenant W. P. Mortenson 
throughout the entire antenna vibration survey; however, separate 
papers are being presented. The purpose of this paper is to 
specifically describe the aeroelastic model used to simulate the 
antenna boom boards, the data handling system utilized in data 
processing for the antenna vibration project, and to present the 
computer program for determining the power spectral density, thereby 
defining the response of the aero elastic wind tunnel model to the 
random aerodynamic forcing function. The detailed analysis and 
explanation of the results provided by this system and recommenda- 
tions for "in the field fixes" for the antenna project can be 
found in [ 1 1 ] . 

The investigation was conducted by the author in conjunction 
with Lieutenant W. P. Mortenson during the period October 1965 
through May 1966 at the U. S. Naval Postgraduate School, Monterey, 
California. My sincere appreciation is extended to the faculty 
and staff of the Aeronautical Engineering Department for their 
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assistance. The author is deeply indebted to Professor L. V. Schmidt 
for his invaluable encouragement and assistance throughout the project. 
2. Description of the Model and Equipment. 

In order to simulate the boom boards in the wind tunnel, it was 
necessary to design an aeroelastic model. The characteristics of 
a dynamic model are quite stringent. Ideally the model should be 
light weight, mounted on frictionless bearings, able to reproduce 
desired frequencies, and sufficiently stiff so as not to introduce 
undesirable structural oscillation. The general approach used was 
to approximate one half of the prototype by hinging the model at 
one end with the other end free, with controllable stiffness. Since 
boom boards are quite long and their end fixity conditions are 
close to simple supports, the first lateral vibration mode was 
assumed dominant. The aeroelastic model can thus be restricted to 
one degree of freedom in the lateral direction. 

Since the intent was to observe many different configurations, 
interchange of the model shape was a necessity. For simplicity, it 
was decided to make the center of the model the strength member. 

An extruded aluminum "I" beam was selected primarily for its light 
weight and ease of machining. This aluminum beam extended the 1 
entire length of the model. 

The outside configuration was formed by shaping balsa wood to 
the desired dimensions. The balsa covering was constructed in two 
symmetric pieces and contoured to fit snuggly on the M I n beam 
(Figure 1). End plates were screwed into the balsa wood at the 
top and bottom; these end plates were mounted to the "I" beam by 
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a bracket. The end plates were adjustable to ensure that the model 
was centered between top and bottom of the tunnel with only a small 
end gap remaining. 

The wind tunnel floor and ceiling were fitted with aluminum 
plates which had the center section recessed slightly to accept 
the end plates. This allowed the model to span the entire height 
of the tunnel and thus reduced the flow around the ends. 

When scaling a dynamic model for wind tunnel testing, the 
variables which must be considered are: fluid density and viscosity, 

fluid velocity, vibration frequency, and inertia effects (stiffness, 
mode shapes, mass distribution, etc,). In low velocity flow the 
fluid may be considered incompressible, i.e., Mach number effects 
may be neglected. 

The first lateral vibration mode has been assumed dominant, 
thus permitting the aeroelastic model to be restricted to one 
degree of freedom in the lateral direction. This allows the 
simplification that model inertia effects need not be scaled with 
respect to prototype. The remaining variables combine into the 
dimensionless parameters: Reynolds number (R=vd/"V ) and 

Strouhal number (S = fd/v). 

Matching of Reynolds numbers between prototype and model 
assures that the aerodynamic forcing function, which causes the 
vibration, correspond. Matching of Reynolds number will effectively 
occur as far as dynamic model response is concerned when the air 
flow about the model is supercritical [ 6 , 7l . Matching of Strouhal 
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number or dimensionless frequency assures dynamic response 
correspondence between prototype and model [8] . 

To illustrate these relations a typical example follows. 



Assume that the prototype bluff body has a characteristic depth 
of 12 inches, a fundamental resonance frequency of five cps, and 



responds to sea level velocities of 25, 50 and 75 miles per hour. 



C!ase Velocity 
mph fps 



S trouhal 
Number 



Reynolds 

Number 



A 25 36.7 

B 50 73.3 

C 75 110.0 



0.136 2.34x10^ 
0.068 4.66x10^ 
0.046 7.00xl(T 



To determine the appropriate aeroelastic model to represent 
the prototype, enter Figure 2 with Reynolds number (R) and model 
characteristic depth (d) to determine the dynamic pressure (q ) . 

With this information enter Figure 3 with q, d, Strouhal number (S) 
to determine the required resonance frequency of the model. Typical 
models which would match the above prototype condition are as follows 



Case 


d 


q 


f 




inches 


psf 


n 

cps 


A 


3- .43 


9.5 


30 




4.0 


14.4 


45 


B 


6.0 


25.2 


20 




5.0 


36.0 


29 




4.0 


57.5 


45 


C 


6.0 


57.5 


20 



Instead of making an exact match on both Reynolds number and 
Strouhal number, certain classes of bluff body cross sections may 
be expected to not be as severely dependent upon Reynolds number 
and hence only require Strouhal number matching. A clue in this 
direction is provided by Figures 5 and 8 of the work by Delany 
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and Sorensen [ 5] where they present steady state drag coefficient 
as a function of Reynolds number for several cross-section shapes. 
The circular cross section is well known to have a transition with 
respect to drag coefficient at a Reynolds number value of about 
3 x 10^ while square and rectangular shapes do not exhibit this 
trait at Reynolds numbers well in excess of 1 x 10 . 

If we presume that the unsteady air loads behave in a *like 
manner to the steady state drag coefficients, then the necessity 
of matching Reynolds number exactly in an aeroelastic model with 
a square cross section loses importance. On this assumption, the 
sample values considered would indicate the necessity to select 
model operating conditions such that Strouhal numbers of 0.03 to 
0.13 corresponded to the model resonance condition for Reynolds 
numbers of approximately 4.0 x 10^. This condition could be 
matched by the model in any of the following configurations: 



d=6 . 0 


inches , 


q= 18.5 


psf , 


f =12 
n 


to 


38 


cps, 


d=5 . 0 


inches , 


q= 21.7 


psf , 


f =16 
n 


to 


48 


cps, 


d=4 .0 


inches , 


q= 42.0 


psf, 


f =28 
n 


to 


80 


cps 



The top of the model, "I" beam, was attached to a rigid support 
by means of a flexure plate (Figure 4). This flexure plate minimized 
frictional damping, while allowing motion primarily in the lateral 
direction Further, it provided some drag restraint. An additional 
drag restraint was provided near the bottom of the model by means 
of a stranded stainless steel cable with turnbuckle arrangement 
for tensioning. 
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The bottom of the model was supported in a manner which provided 
the c&f^>iHt$ of varying the models natural frequency. A heat 
treated one-half inch diameter steel drill rod was attached to the 
aluminum "I" beam by two bolts machined to accept ball bearings. 

The bolts were adjustable to allow centering of the drill rod, and 
the ball bearings provided a relatively frictionless pin joint 
(Figure 5). The drill rod was cantilevered to a movable steel 
support block, precision drilled to accept the spring rod, and 
provided with six set screws which fixed the bar in position. The 
movable steel support block was fitted to a 16 inch precision 
machined wedge shaped track. Eight lag bolts secured the sliding 
block in any position (Figures 6 A,B,C). The track was bolted to 
a one-inch thick steel plate supported by four welded four-inch 
diameter steel pipes which formed legs in a manner similar to a 
dynafocal suspension principle (Figure 7). This rebuilt model 
support base increased the rigidity of the model's reference base 
over that used in [10]. The improvement in the lateral stiffness 
constraint provided a means for more rapid adjustment of the spring 
bar length while assuring that the base of the constraint was attached 
to the support base. Figure 8 illustrates the complete specifications 
of the model. 

The wind tunnel utilized was the 42 inch by 60 inch subsonic 
tunnel located at the U. S. Naval Postgraduate School, Monterey, 
California. This tunnel is a closed section, vertical return, 

200 mile per hour installation built by West-Coast Research 
Corporation of Los Angeles, California. Power is supplied by two 
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150 horsepower electric motors driving two independent constant 
speed variable pitch fans. However, for ease of operation only one 
fan was utilized; this fan provided adequate power for the dynamic 
pressure required. In addition, a survey of the dynamic pressure 
on tunnel centerline was completed. In order to check the turbulence 
level of the tunnel, a hot-wire survey was performed. The results of 
these surveys are presented in Appendix A. The aeroelastic model 
was mounted vertically in the test section as shown in Figure 9 
along with an illustration of the entire data handling system. 

The primary instrument for response measurement was the Bently 
displacement measuring device. This device provides an instanta- 
neous measurement of the lateral displacement of the bottom of 
the beam. 

The Bently sensor unit, a small pancake shaped coil, is 
excited at a one megacycle carrier frequency by a transistorized 
amplifier unit and power supply combination. The rapidly 
fluctuating magnetic field at the face of the coil allows sensing 
of the air gap change between the coil face and a flat non-magnetic 
conducting surface due to eddy current effect. The frequency limit 
of the Bently is 10,000 cycles per second. This is well above the 
125 cycles per second anticipated during the antenna survey. In 
practice this device proved quite reliable with a consistent distance 
resolution of about 20 micro-inches. Figure 6B shows the Bently 
sensor unit mounted to the model support base. Therefore, the gage 
measured relative displacement between the support floor and the 
model. Appendix B provides a calibration for the Bently inductance 
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displacement measuring device. For test purposes, the air gap was 
allowed to vary in the linear range from 0.02 to 0.08 inches, where 
the output sensitivity was approximately 70 volts/inch. 

The secondary source of model frequency response data was the 
Gulton type A321 accelerometer. This piezoelectric type device was 
used with an Endevco amplifier. The accelerometer has a flat 
response curve to above 3,000 cycles per second. Results of Gulton 
accelerometer calibrations are presented in Appendix B. It was 
decided to use the Bently as the primary source of data due to its 
reliability and ability to measure relative displacement between 
fixed reference and the base of the model instead of the inertial 
reference frame of the accelerometer (Figure 6C). 

Figure 10 presents a block diagram of the entire data 
acquisition system. Due to the complexity of the wiring and the 
desire to have an orderly system, BNC fittings were used throughout 
in conjunction with shielded coaxial cable. The coaxial cable 
decreased the possibility of picking up any outside disturbances 
such as 60 cycle from power supplies, ground loops, and 
f lourescent lights . 

All data were monitored and controlled through a central switch 
panel in a manner similar to a telephone patch panel in order to 
maintain a systematic control of data flow. Voltage signals from 
the Bently and accelerometer were monitored and recorded using a 
Ballantine model 320A voltmeter in order to obtain the true RMS 
reading of the randomly varying signal. This meter uses a diode 
squaring circuit so as to give an accurate measure of the RMS 
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voltage, and can handle situations where peak values occur up to 
ten times greater than the RMS . 

The Vidar, integrating digital voltmeter, reads both frequency 
and voltage. It was utilized in the system to insure that the 
initial setting on the Bently, with the model stationary, was in 
the center of the linear range of the response characteristics of 
the Bently. Its secondary function was to calibrate the center 
frequency of the tape recorders frequency modulating (FM) 
data channels. 

The tape recorder /reproduce unit used was the Ampex FR-100. 

This provided 14 channels for the recording of data. Both frequency 
modulated and direct record amplifiers were available to record 
data. With a recorder /reproduce speed of seven and one-half inches 
per second, the frequency response for the frequency modulated 
amplifier was from zero to 1250 cycles per second. The input to the 
recorder was limited to one volt RMS or * 1.40 volt peak to peak. 

The voltage signals from the Bently and accelerometer were recorded 
through the frequency modulated amplifier. Voice announcements 
were recorded through a direct response amplifier. These 
announcements properly labeled and identified the data. The system 
recorded the data on one inch mylar instrument grade tape so that 
the voltage signals which represent an analog of the model response 
could be analyzed at a later date. 

The ifcektronic 545A cathode ray oscilloscope is a single beam 
device which behaves in a dual beam manner using the dual trace 
amplifier channel, and it was used to monitor various signals. 



22 



The dual trace provided the means for simultaneously observing the 
signals from both the Bently and the accelerometer. Since the input 
to the Ampex tape recorder was limited to one volt RMS, observing 
the input determined whether attenuation was necessary to keep the 
signal within limits. Finally, the scope could be connected to the 
Ampex reproduce amplifier (output) to insure that the data was being 
properly recorded and reproduced. 

3. Test Procedures. 

The entire data acquisition system with a few exceptions 
consisted of tube elements. Therefore, prior to obtaining data, the 
entire system was allowed a proper warmup time in order to insure 
that all thermal transients had settled out. The Ampex tape 
recorder/reproduce amplifiers were calibrated in accordance with 
the operation manual. 

Calibrations were recorded on the magnetic tape through the 
frequency modulated amplifiers on the exact channels that the 
Bently and acceleration voltage signals were to be recorded. This 
consisted of zero volts and plus and minus 1.4 volts (DC). This 
calibration could be utilized later during the data analysis phase. 

To determine the system damping and to calibrate the spring rod 
length as a function of the natural frequency of the model, static 
tests were performed with the wind tunnel inoperative. The model 
base was mechanically displaced from the centered equilibrium 
position and sharply released. The model’s decaying response of 
both displacement and acceleration was recorded onto magnetic tape. 
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The vibrational response of the model was obtained by a series 
of test runs. Each was run at a fixed natural frequency set by 
positioning the calibrated spring rod to the proper length while 
varying the tunnel dynamic pressure. Here again, the data from 
the Bently and accelerometer were recorded on one-inch magnetic 
tape and the RMS values from the Ballantine voltmeter were read and 
recorded in the data logs. Visual observations, using the Tektronic 
oscilloscope, were made of the response wave form to insure that 
the amplitude did^not exceed record amplifier limitations, and 
further to check that data was being accurately recorded 
and reproduced. 

Voice announcements on a separate channel preceded and followed 
the calibration check, model frequency response check and the 
dynamic data runs. This insured that the data was properly labeled 
and would definitely be retrievable. 

The data from the one -inch magnetic tape was next played back 
through an analog representation of an averaging circuit and plotted 
to observe the continuity of the curves, and further to look for any 
violent or sporatic portions. As expected from the oscilloscope 
observations, the curves were smooth and well-behaved, although 
complex and random in nature. Further, it appeared that the data 
was a stationery random process (this was subsequently confirmed 
by power spectral analysis.) Since the data was well-behaved, any 
portion of the data could be analyzed without fear of large errors. 
The model frequency response check and portions of the data were 
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plotted on strip chart Brush oscillograph to obtain frequency. 
Samples of this data may be found in Appendix C. 

4. Data Processing System. 

At this point, approximately 450 data runs in analog form 
were recorded on one-inch magnetic tape. It was desirable to form 
the auto-correlation function and then the power spectral density 
in order to obtain the response traits of the various configurations 
to the random forcing function. It was estimated that it would take 
approximately three to five man-months to process each run by hand. 
Obviously it was necessary to utilize the Control Data Corporation 
(C.D.C.) 1604 digital computer. However, a major hurdle still 
existed, namely the conversion of the analog data into digital form, 
acceptable as computer input. 

The process of converting real analog data into computer 
digital form was accomplished in two distinct and separate steps. 

In the first step the analog data was digitized at discrete points 
and written in BIT form on one-half inch magnetic tape. The second 
step was to unpack the digital data and put it into a form useable 
in the C.D.C. 1604 as standard program input. 

Figure D-l is a block diagram of the analog to digital 
conversion system. The wind tunnel data from the one -inch magnetic 
tape was reproduced through the Electrical Engineering Departments 
Ampex CP-100, a solid state record/reproduce system. The digitizing 
computer was the CDC-160 located in the Electrical Engineering 
Department Computer Laboratory. The CDC-160 selects either the 
12 BIT analog to digital converter unit for input, or the CDC-163 
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tape unit for output. The computer has a 4096 cell 12 BIT memory 
length. The lower 4,000 cells were reserved for storage, and the 
top 96 cells were reserved for programming. The core storage was 
an obvious limiting factor. However, 4,000 data samples were 
assumed adequate. The basic concept of the operation was to take 
samples until the core storage was full (4,000 samples) and then 
write this information onto one-half inch magnetic tape. 

The procedure utilized a manually initiated pulse to trigger 
the analog to digital converter. The computer controls the inter- 
sample time delay and the analog to digital converter samples rapidly 
until 4,000 discrete samples have been taken and stored in the memory 
core. The computer then dumps the output onto one-half inch magnetic 
tape. The analog to digital converter then remains hung up until 
receipt of another manual pulse. Program Digitize, a CDC-160 tape 
program, is found in Table D-l. 

A sample rate of 500 discrete points per second allows spectral 
distributions up to 250 cycles per second, while a digital record of 
eight seconds duration provides a lower bound on spectra of 
approximately one cycle per second. This frequency range is more 
than adequate since the natural frequencies of aeroelastic models 
tested varied from 35 to 90 cycles per second. 

In order to identify each block of 4,000 data points on the 
tape, a 16 digit identifier was first written followed by the data 
points. These identifiers are extremely important since they are 
the only means by which the data may be retrieved at a later time. 
This 16 digit identifier is a manual input to Program Digitize. 
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The data in digital form may be recalled by subroutine 
Data. Subroutine Data, Table E-2, in turn uses subroutine Findit, 

Table E-3, which locates the identifier on the tape and subroutine 
Unpack, Table E-4, a machine language routine to unpack the samples 
so that one-to-one correlation exists between the original CDC-160 
listing and the unpacked CDC 1604 version. 

A short CDC-160 tape program Test, Table D-2, was used to 
adjust the bias voltage to minus five volts and the desired inter- 
sample delay to the proper time. Additional more detailed 
information of the digitizing system can be found in Lieutenant 
N. S. Barrett 1 s thesis "Extracting Analog Signals From Noise Using 
The Digital Computer" [ 9] • 

In order to test the capabilities of the digitizing system, a 
test tape was prepared with a 250 cycle per second sine wave, square 
wave, triangular wave, 20 cycle per second band-limited white noise 
and finally white noise. The data was analyzed and presented 
graphically, Figures D2-6. The system worked successfully and 
provided very satisfactory results. 

The final step in the data handling system was the computer 
program, written in Fortran 63, which utilized the Tukey method 
for estimating the power spectral density [4] . 

The complete program may be found in Appendix E. The program* 
possesses sufficient comments to make it self-explanatory. However, 
some additional comments on the input cards are provided in Figure E-l. 
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Basically, this program used the data recalling system 
previously described. The integer format^ from the subroutines 

r- 

was ^converted to floating point format and converted into volts 
by a constant. In order to insure that no bias existed, the data 
were normalized to wipe out the DC term on an average basis. 

The auto-correlation function of the random response was 
computed, and subsequently used, to estimate the spectral distribution 
of the energy by applying a Fourier cosine transform to the auto- 
correlation function. For any given signal, aperiodic or random in 
nature, it is quite clear that the auto-correlation function and 
spectral density are quite meaningful, but only the auto-correlation 
function is directly measurable. The power spectral density may be 
defined as the energy distribution on an average basis in the 
frequency domain; that is, energy with respect to frequency. 

Appendix F gives a more detailed mathematical explanation of the 
auto-correlation function and power spectral density to show the 
inter-relationships. The final step in the program was to determine 
the energy contained in the power spectrum by application of the 
trapezoidal rule. Comparing this result with the total energy 
provides a rapid indication of the accuracy of the results obtained. 

A sample listing of the data output is given in Table E-5 and is 
graphically presented in Figure E-2. The data and graphs specifically 
define the frequency response of the aeroelastic model to the 
aerodynamic random forcing function. This information is imperative 
to the investigation of the aerodynamic phenomena and subsequent 
analysis described in [ill. 
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5. Summary. 

A controlled test environment for conducting wind-tunnel 
investigations of bluff body aeroelastic models has been established. 
The data acquisition and processing systems are precise, orderly, 
and operational. The conversion from an analog voltage record to a 
digital format has been accomplished with a high degree of accuracy. 
The computer program was utilized to analyze approximately 450 wind 
tunnel data runs, and presented the output in numeric and graphic 
form. These results were imperative to the success of the 
investigation and analysis of aeroelastic phenomena conducted on 
the various configurations as presented in [ 11] . 
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FIGURE 1 
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FIGURE- 4 
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FIGURE 6 A 
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FIGURE 6B 
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FIGURE 6 C, 
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APPENDIX A 



WIND TUNNEL TEST SECTION CALIBRATION 



The test section pressure survey shown in Figure A-l was 
performed at the U. S. Naval Postgraduate School's low speed wind 
tunnel. The profiles shown were taken 22 inches aft of the test 
section leading edge, with reference dynamic pressure of 40 pounds 
per square foot - 

These wind tunnel flow traits were improved with respect to 
tunnel dynamic pressure and flow inclination. Although the 
importance of these changes upon bluff body separation may not be 
critical or discernible, it is good practice to have a well 
established and uniform test environment in the wind tunnel 
test section. 

Aconstant temperature hot wire was placed on tunnel centerline 
at the model test station in order to make clear tunnel turbulence 
measurements. The hot wire was 0.001 inch diameter tungsten, 
approximately 0.08 inches in length and was supported on the end 
of two needles projecting from a support tube. It was possible 
to make direct measurements of turbulence level relative to free 
stream velocity using a prototype constant temperature amplifier 
circuit recently developed at the U. S. Naval Postgraduate School. 
These turbulence measurements were read directly by the Ballantine 
True RMS voltmeter, and also recorded on magnetic tape for 
subsequent digital analysis using procedures described elsewhere. 
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The direct measurements of turbulence indicated a turbulence 



level of approximately 0.45%, but the single wire sensed 
fluctuations both in the vertical and axial direction. On the 
assumption of isotropic turbulence, it is possible to estimate 
the axial velocity turbulence level as approximately 0.32%. This 
value was obtained uniformly for a test section q range of 10 to 
40 psf. The above value of turbulence level is considered high 
for high calibre wind tunnel work in a low turbulence environment, 
but considering that no anti -turbulence screens were installed 
in the settling section, and that the test was concerned with 
bluff body separation, the value was quite acceptable. The 
installation of screens in the settling section should result in a 
much lower value of turbulence level. 

The results of the digital calculations are shown on Figure A-2 
with respect to both auto-correlation and power spectral density. 

The area of the mormalized auto-correlation curve is 0.102 seconds 
which may be interpreted as a measure of the scale of turbulence, 
that is, on the average, turbulent fluctuations are unrelated 
for time intervals beyond 0.102 seconds. The power spectral 
density curve shows no dominant energy peaks other than one at 
1/2 cps . It is quite possible that the 1/2 cps peak was caused 
by the vibration of the hot wire support. The absence of other 
peaks indicates an absence of noise due to fan blade rotation or 
an undamped flow separation from other sources in the tunnel 
circuit . 
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The basic conclusion reached from the hot wire measurements 
is that the test section circuit is functional in its present 
configuration and the addition of screens at a later date will 
make the facility very useful for fundamental low turbulence 
experiments . 

The tunnel airspeed and/or dynamic pressure was set using 
clear tunnel values. In making a final evaluation, ‘solid blockage 
effects may be taken into account. An estimate of blockage for 
the case of a cylinder spanning the tunnel is given by the formula 

^ - 7 ? a. 

£ S TJX 



where a = radius of the cylinder 
h = height 

c = fractional velocity increase 

O 



model width 
inches 

4 

5.7 

8 



.0146 

.0296 

.0584 



'corr^ref 

1.03 

1.06 

1.12 
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FLOW SURVEY 

3.5 X 5,0 FT. L.OW SPEED WIND TUNNEL 
- CLEAR TEST SECTION 
. q = 40 PS F AT TEST SECTION l.E. 

t BE F, 



CEILING CEILING 





FIGURE A-l 

Surnmqry of Wind Tun nef Flow Characteristics 
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APPENDIX B 



CALIBRATION OF BENTLY AND GULTON INSTRUMENTS 

The Bently sensor head was mounted in an aluminum test 
block. The procedure used was to obtain the gage voltage for 
a specific distance measured by a micrometer. The plot of 
distances versus output voltage is shown in Figure B-l. The 
slope of the linear region was 70 volts/inch. The Bently sensor 
unit was serial number 14184, and the distance detector unit 
model D-252 was serial number 5254. The ease of calibration and 
the consistent reliability made the Bently displacement measuring 
system an extremely useful instrument. 

The Gulton model A321, a piezoelectric accelerometer, is 
extremely sensitive and because of its high natural frequency 
it can be used through a wide frequency range. However, since 
it is an A.C. instrument it cannot be statically calibrated. 
Therefore, the accelerometer was calibrated against a previously 
calibrated Statham on the Calidyne shaker table. The Gulton 
accelerometer serial number 1555 was combined with an Endevco 
amplifier serial number 128. The calibration is shown on 
Figure B-2. 
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FIGURE B-2 



CALIBRATION OF GULTON ACCELEROMETER 
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APPENDIX C 



BRUSH OSCILLOGRAPH TRACES OF DATA 
AND FREQUENCY RESPONSES 

Typical brush recorder tracings of response data of the 
aeroelastic model to the random forcing function are shown on 
Figure C-l. The lever arm and average frequency are shown on 
each trace. 

Typical brush recorder tracings of the static frequency 
response test are shown in Figure C-2. These are used to determine 
the system natural frequency (f^) and damping ( p ) . 



48 



FIGURE C-l 



v 

TYPICAL BRUSH RECORDER TRACINGS OF RESPONSE DATA 
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FIGURE C-2 



TYPICAL BRUSH RECORDER TRACINGS OF 
STATIC FREQUENCY RESPONSES 



L = 9.5" f = 20.75 cps /3 = 0.0175 

b n 





1, = 6.0" f = 34.5 cps /3 Poorly defined 

b n 
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APPENDIX D 



ANALOG TO DIGITAL CONVERSION SYSTEM 



The process of converting real analog data into computer 
digital form was accomplished by the system shown in Figure D-l. 
The system utilized the C.D.C. 160 computer with Program 
Digitize, Table D-l, to sample the analog data digitize and 
write the digital form on magnetic tape compatible with the 
C.D.C. 1604 computer. Program Test, Table D-2, is used to 
calibrate the system to insure accuracy. The data in digital 
form may be recalled by utilizing subroutine Data, Table E-2, 
which in turn utilizes subroutine Findit, Table E-3, which 
locates the pre-set identifier on the tape and subroutine 
Unpack, Table E-4, a machine language routine which unpacks 
the data. In order to test the capabilities of the digitizing 
system a test tape was prepared The first 160 points of each 
is plotted in Figures D-2A through D-2E. 
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TABLE D-l 



PROGRAM DIGITIZE 



0 


0101 


PTA 




1 


0603 


ADN03 


Jump ^to 


2 


7064 


JPI64 


Initial 


3 


7500 


EXFOO 




4 


2410 


2410 


Set enable 


5 


7500 


EXFOO 




6 


1401 


1401 


A/D Ch. 1 


7 


7600 


INA 




10 


4176 


STI76 




11 


2076 


LDD76 




12 


3465 


SBD65 




13 


6134 


NZF34 


60 to 47 


14 


7500 


EXFOO 




15 


2400 


2400 


Clear enable 


16 


2074 


LDD74 ) 




17 


4160 


STI60 ) 




20 


2067 


LDD67 ) ) 


Form indent 


21 


4161 


STI61 ) 




22 


2075 


LDD75 ) 




23 


4162 


STI62 ) 




24 


2073 


LDD73 ) 




25 


4163 


STI63 




26 


7500 


EXFOO 




27 


2111 


2111 


163 T. V. 


30 


7303 


OUT03 


Output From 


31 


0000 


C 0000 


LWA + 1 


32 


6102 


NZF02 




33 


0133 


0133 




34 


2246 


LDF46 


0137 A 


35 


4076 


STD76 




36 


2075 


LDD75 




37 


3473 


SBD73 




40 


6155 


NZF55 


G0 to 0115 


41 


0401 


LDN01 




42 


4075 


STD75 


Set block No. = 1 


43 


5455 


A0D55 




44 


3454 


SBD54 




45 


6153 


NZF53 


G0 T0 0120 


46 


6061 


ZJF61 


G0 TO D127 


47 


5476 


A0D76 ) 




50 

51 


2066 

0701 


LDD66 ) ) 

SBN01 ) 


Inter sample 
delay loop 


52 


6501 


NZB01 ) 




53 


7056 


JPI56 ) 




54 


0760 


0760 


Max blocks 760 


55 


0000 


0000 


Total block count 
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56 

57 

60 

61 

62 

63 

64 

65 

66 

67 

70 

71 

72 

73 

74 

75 

76 

77 

100 

101 

102 

103 

104 

105 

106 

107 

110 

111 

112 

113 

114 

115 

116 

117 

120 

121 

122 

123 

124 

125 

126 

12 7 

130 

131 

132 



TABLE D-l (Con't. ) 



0005 


0005 




0003 


0003 




0133 


0133 


1st word address 


0134 


0134 


2nd of header 


0135 


0135 


3rd of address 


0136 




LWA of header 


0100 




Addr of initial 


0000 


C 


LWA of data 


0000 


M 


Inter sample delay 


0000 


M 


Spare ID (anything) 


0000 


M 


Initial run no 


0000 


M 


No runs wanted 


0000 


M 


No samples/block 


0000 


M 


No blocks/run 


0000 


C 


Current run number 


0000 


C 


Current block number 


0000 


c 


Running Storage 
Address 


0000 


c 


Jump control 


4077 


STD 7 7 


Enter initial 


2200 


LDC00 




0137 


0137 




4076 


STD 76 


Set FWA of data 


3072 


ADD 72 


Compute LWA + 1 


4031 


STD31 




0701 


SBN01 




4065 


STD65 


Set LWA of data 


2070 


LDD70 


Load initial run no. 


4074 


STD 74 


Set I. R. N. 


0401 


LDN01 




4075 


STD75 


Set I. B. N. 


7077 


JPI77 


Exit initial 


5475 


A0D75 


Update block no. 


5455 


A0D55 


Update block count 


7057 


JPI57 


Go to 0003 


2074 


LDD74 




3471 


SBD71 


Check no . runs 
enough 


6103 


NZF03 




2074 


LDD74 


Current run no. 


7701 


MLT01 




5474 


A0D74 


Update run no. 


7057 


JPI57 


Go to 0003 


7500 


EXFOO 




1111 


1111 




2055 


LDD55 


Display total 
block count 


7702 


HLT02 
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TABLE D-l (Con't . ) 



133 




0000 


Run no. 10 


134 




0000 


Spare ID 


135 




0000 


Block ID 


136 




0000 


No. blocks/run ID 


137 


) 








\ 


) Data storage area 




7776 


) 







Note: M = Manual entry 

C = Computed entry 
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0 

1 

2 

3 

4 

5 

6 

7 

10 

11 

12 

13 

14 

15 

16 

17 

20 

21 

22 

23 

24 

25 

26 

27 

30 

31 

32 

33 

34 

35 

36 

37 

40 

41 

42 

43 

44 

45 



TABLE D-2 



PROGRAM TEST -160 



7500 


EXFOO 


2410 


2410 


7500 


EXFOO 


1401 


1401 


7600 


INA 


4070 


STD70 


7500 


EXFOO 


2401 


2401 


7303 


OUT03 


0072 


0072 


6102 


NZF02 


0070 


0070 


0400 


LDN00 


4071 


STD71 


6414 


EJB14 


7700 


HLTOO 


7500 


EXFOO 


2410 


2410 


7500 


EXFOO 


1401 


1401 


7600 


INA 


4143 


STI43 


0300 


N0P 


0300 


N0P 


0300 


N0P 


0300 


N0P 


0300 


N0P 


0300 


N0P 


0300 


N0P 


0300 


N0P 


2042 


LDD42 


0701 


SBN01 


6501 


NZB01 


7044 


JPI44 


M 


Timing 


0043 


0043 


0022 




7701 


HLT01 



Set enable 
Call A/D CH 1 

Output 

From 0071 to 0072 

Return to 000Z 
Set enable 
Call A/D CH 1 



control word 
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FIGURE D-l 



DIAGRAM OF ANALOG TO DIGITAL CONVERSION SYSTEM 
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FROM DIGITAL TEST TAPE 



FIGURE D-2B 
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250 CYCLE PER SECOND TRIANGULAR WAVE 
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FROM DIGITAL TEST TAPE 
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FROM DIGITAL TEST TAPE 
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FROM DIGITAL TEST TAPE 



APPENDIX E 



POWER SPECTRAL DENSITY PROGRAM 



The energy distribution on an average basis in the 
frequency domain is defined by the power spectral density. A 
computer program utilizing the Tukey method for estimating 
spectral distribution was written - Table E-l. Tables E-2 
through 4 are subroutines which locate and unpack the data 
desired. The program has sufficient comments to make it self- 
explanatory. Figure E-l is an illustration of the control card 
utilized in the program. A sample listing of the data output is 
given in Table E-5 and graphically presented in Figure E-2. The 
data and graphs specifically define the frequency response of the 
aeroelastic model to the random aerodynamic forcing function. 

It will be noted in Table E-5 that the first 100 values of 

input data are listed in order to give visual assurance that the 

numbers are valid. The listing of auto -correlation is normalized 

with respect to its value at Z = 0 and is specified as AN(TAU) 

versus Tau(SEC). The value A(0) shown above the listing is the 

actual mean square value of the auto -correlation at time shift 

2 

/ C*= 0 seconds, and is in units of (volts) . This value may be 
related to the tape recorder voltage signal by the scaling factor 
of 5.0 volts (digitized) = 1.4 volts (recorded). The power 
spectral density is listed as XN(FREQ) versus FREQ(CPS) and is 
normalized such that the area under the curve of XN versus FREQ 
is unity. 
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TABLE E-l 

PO.VER SPECTRAL DENSITY COMPUTER PROG RAJ/ 

-COOP, ,JESBERG BOX J* I /I /0/49/S/l S/2S/E/45*54/ 5= 50/6*51 ,15, 20000* 
-FTN*L*E. 

PROGRAM POWSPEC 

C TUKEY METHOD FOR ESTIMATING POWER SPECTRAL DENSITIES 

DIMENSION P ( 4000 )*A(501),X(50l) *FREQ( 501 ) *TAU( 501 ) *Q( 3000) 
DIMENSION KDATA(4000,2 ) , I DENT ( 200 ) 

DIMENSION I T I TLE( 12 ) 

DIMENSION DATE( 160),SKIP( 160) 

COMMON KDATA 
EQUIVALENCE (KDATA, P) 

FIRST STEP IS DATA CONTROL CARD— DATA READ IN BY TAPE OR CARD 
FIND IF DATA IS ON TAPE — ITAPE=0 FOR TAPE*=1 FOR CARD DATA 
WRITE (6,103) 

1 READ (5,101) I TAPE , IDELTA »DELT AT»ISTART» I MAX ,M 
WRITE (6,104) 

WRITE( 6,101) ITAPE.IDELT A *DELT AT , I START , I MAX ,M 
IQ=I START-1 
I ADJ= I MAX-IQ 

2 I F ( ITAPE - 1) 6,5,5 

START DATA CARD READ IN — USE THIS BRANCH FOR TEST PURPOSES 
IF CARD DAJA REMOVE STATEMENT 11 AND REPLACE WITH A CONTINUE 

5 READ (5,102) ( P ( I ) , I =1 , IMAX ) 

WRITE (6,105) 

WRITE (6*106) (P( I ), 1=1*100) 

GO TO 30 

START TAPE DATA READ IN — USE THIS BRANCH FOR PRODUCTION 
NUMHDR * NO. OF HEADINGS 

6 READ 1 10.NUMHDR 
DO 7 N=1 *NUMHDR 

7 READ 111 * IDENT ( N) 

WRITE (6*114) 

WRITE (6,111) ( IDENT(N) ,N=1, NUMHDR) 

C MAIN DO LOOP ON ENTIRE PROGRAM 

DO 500 N=l, NUMHDR 
IDENT =IDENT (N) 

KL I ST= 1 
MAX= IMAX 

CALL DATA( IDENT, MAX, KLIST»KFLAG) 

IF(KFLAG) 8,9*8 

8 WRITE (6*112) 

WRITE (6.113) N* IDENT 
GO TO 500 

9 DO 11 J*1 * IMAX 
1 = 1 

10 P( J ) *KDATA( J* I ) 

11 P(J)=P( J)/(-409.6) 

WRITE (6*105) 

WRITE (6*106) (P( I )* 1-1*100) 

60 TO 30 
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TABLE E-l (Con't.) 



i 



C NORMALIZE INPUT DATA TO WIPE OUT D-C TERM ON AVERAGE BASIS 

30 ASUM=0 • 0 

31 DO 32 I - l ♦ I ADJ * IDELT A 
ASUM=ASUM+P( I ) 

32 CONTINUE 

C AZ IS NUMBER OF DATA SAMPLES BEING COMPUTED 

AZ-FLOATF ( ( IADJ-l)/( I DELTA )) +1.0 
ASUM- ( ASUM) / { AZ ) 

33 DO 34 I =1 » I ADJ * IDELTA 
P( I ) =P ( I ) -ASUM 

34 CONTINUE 
ASIS-0.0 

DO 35 1=1 » 160 
DATE ( I )=P( I ) 

SKIP! I )*ASIS 
AS I S-AS IS + 1 *0 

35 CONTINUE 

C FIND MEAN SQUARE AZERO I.E. AUTOCORRELATION AT TAU-ZERO 

40 ASUM-0 *0 

41 DO 42 I=1*IADJ*IDELTA 
ASUM-ASUM+PI 1 )**2 

42 CONTINUE 
AZERO-ASUM/AZ 

C FIND AUTO-CORRELATIONS A( L ) FOR L«1*M 

43 DO 46 L = ltM 
ASUM-0 . 0 
BSUM-0.0 
CSUM-0.0 

MZ* ( L* I DELTA+1 ) 

44 DO 45 I -MZ * I ADJ * IDELTA 
I Z= I -( L*IDELTA) 

ASUM=ASUM+P( I Z ) *P ( I ) 

BSUM=BSUM+P( IZ) 

CSUM-CSUM+P ( I ) 

45 CONTINUE 

AZ*( 1.0)/FLOATF( ( ( IADJ-1) /( IDELTA) )+l-L) 

A( L ) -( AZ ) *( ASUM ) - ( AZ**2 ) *BSUM*CSUM 

46 CONTINUE 

C FIND XZERO POWER SPECTRAL DENSITY AT FREQ-ZERO 

50 ASUM-O* 0 
MZ-M-1 

FM-FLOATF ( M ) 

CSl-COSFI 3*14159/FM) 

SN1-SINFI3.14159/FM) 

CSL-CS1 

SNL-SN1 
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TABLE E-l (Con't.) 



51 DO 52 L=1 »MZ 
AZ=HANNING FACTOR 

FIND COSINE SUM FROM TRIG# IDENTITIES 
AZ*1 #0+CSL 
ASUM=ASUM+AZ*A(L» 

CSL1=CSL*CS1-SNL*SN1 
SNL1=SNL*CS1+CSL*SN1 
CSL=CSL1 
SNL=SNL 1 

52 CONTINUE 
DZ*0.5/FM 

XZERO=DZ*( ASUM+AZERO) 

C FIND X ( K ) POWER SPECTRUM AT K«1»M 

CSK=CS1 

SNK=SN1 

53 DO 59 K*1.M 
ASUM«0,0 
CSKL*CSK 
SNKL*SNK 
CSL*CS1 
SNL=SN1 

54 DO 55 L*1»MZ 

AZ = ( 1 • 0+CSL ) *CSKL 
ASUM=ASUM+AZ*A(L) 

CSL1=CSL*CS1-SNL*SN1 
SNL1=SNL*CS1+CSL*SN1 
CSL=CSL 1 
SNL*SNL 1 

CSKL1*CSKL*CSK-SNKL*SNK 

SNKL1=SNKL*CSK+CSKL*SNIC 

CSKL=CSKL1 

SNKL=SNKL1 

55 CONTINUE 

I F ( K-M ) 56 *57 *57 

56 DZ= 1 #0 
GO TO 58 

57 DZ=0#5 
GO TO 58 

58 DZ=DZ/FM 

X ( K ) =DZ* ( ASUM+AZERO ) 
CSK1=CSK*CS1-SNK*SN1 
$NK1*SNK*CS1+CSK»SN1 

csk=cski 

SNK-SNK1 

59 CONTINUE 
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TABLE E-l (Con’t.) 



SHIFT INDICES— ZERO FREQ WILL START AT K»1 AFTER SHIFT 
USE STORAGE SPACE OF P(I) DATA DURING SHIFT OPERATION 
NORMALIZE ANUTOCORRELAT ION 
NORMALIZE AUTOCORRELAT IONW.R*T* AZERO 

60 DO 61 1*1, M 
P( I+1)=X( I ) 

P< 1+1001 ) = A( I ) 

61 CONTINUE 
X< 1 )=XZERO 
A( 1 )*1.0 
MZ=M+1 

62 DO 63 1=2, MZ 
X( I)=P( I) 

A( I )=P( 1+1000 ) / ( AZERO ) 

63 CONTINUE 

APPLY TRAP* RULE TO FIND ENERGY CONTAINED IN POWER SPECTRUM 

FOR RANGE OF K=1»M+1- DEFINE ENERGY AS XENGY 

NOTE X C K ) IS ENERGY (W*R.T* UNIT CHANGE OF INDEX K 
ASUM=0 • 0 

65 DO 66 K=2*M 
ASUM= ASUM+X ( K ) 

66 CONTINUE 

XENGY=0* 5*( X ( 1 ) +2* * ASUM+X ( M+l ) ) 

C FIND FRACTION OF TOTAL ENERGY IN CALCULATED FREQ* RANGE 

X FACT 3 XENGY /AZERO 

C OBTAIN SPECTRAL DENSITY W*R*T* FREO(CPS) — NORMALIZE W*R*T XENGY 

FREQ ( 1 ) =0 *0 
T AU ( 1 ) =0* 0 
MZ=M+1 

67 DO 68 K=1,MZ 

AZ= 2.0*(DELTAT)*FM 

X ( K ) 3 ( AZ ) *X ( K ) /XENGY 

FREQ ( K+l ) =FREQ( K) + ( 1 . 0 ) / ( AZ ) 

TAUIK+l ) = TAU( K) +DELT AT 
60 CONTINUE 



\ 
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C POWER SPECTRUM WRITE OUT INSTRUCTIONS 

70 WRITE (6,103) 

WRITE (6,107) XFACT , AZERO 
WRITE (6,108) 

WRITE (6,109) (TAU(K) *A(K) »FREQ(K) *X(K) ,K»1»MZ) 

LABEL=4H 

ITITLE( 1 ) *8HJESBERG 
I TITLE ( 2)=8HPOWER SP 
ITITLE( 3)=8HECTRAL D 
ITITLE(4)=8HENSITIES 
ITITLE(5)«8H — RUN 
READ (5,115) I T I TLE( 6 ) 

DO 80 J=7 » 12 
80 ITITLE( J) =8H 

CALL DRAW ( 125 »FREQ »X ,0 ,0 , LABEL , I T I TLE* 20 • , .2 ,0 *0* 2, 2 ,8 ,2 , 1 »LAST ) 
500 CONTINUE 

101 FORMAT (212, F8#5, 216,13) 

102 FORMAT ( 12F6*4) 

103 FORMAT (99H1 TUKEY ^SPECTRUM ESTIMATES* 

1 POWER SPECTRAL DENSTIY AND AUTO-CORR# CALCS# // ) 

104 FORMAT (20H DATA CONTROL CARD //) 

105 FORMAT (30H0 DATA RECORD — (P(I)«1*100) //) 

106 FORMAT (10F8#3) 

107 FORMAT (8H XFACT= F8#5,8H A(0)« F12#5 //) 

108 FORMAT (42H TAU(SEC) AN(TAU) FREQ(CPS) XN(FREO) // ) 

109 FORMAT ( F9# 3 ,F 10# 5 , F12 #3 ,F10#5 ) 

110 FORMAT ( 3X , 1 4 ) 

111 FORMAT ( 3X» 016 ) 

112 FORMAT ( 31H1 ERROR IN SR DATA AT HEADING ) 

113 FORMAT (40X,I3,7X,O16) 

114 FORMAT ( 22H IDENT ( N ) *N«1 »NUMHDR // ) 

115 FORMAT ( 3XA8 ) 

END 
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TABLE E-2 



SUBROUTINE DATA ( 1 DENT »MAX ♦ KL I ST » KFLAG ) 

DIMENSION I BLOCK. (1001 ) , KDATA < 4000 » 2 ) 

COMMON KDATA » I BLOCK 

CALLING INSTRUCTIONS FOR SUBROUTINE DATA 

ARG 1. IDENT IS A 16 OCTAL DIGIT IDENTIFYING NUMBER 
SUPPLIED BY THE 160 PROGRAM WHEN TAPE WAS MADE* 

THE CALLING PROGRAM MUST SUPPLY THIS NUMBER IN ORDER FOR 
THE CORRECT BLOCK OF DATA TO BE RECALLED. 

ARG 2. MAX IS THE NUMBER OF SAMPLES PER BLOCK AND MUST BE LESS 
THAN 4000 DEC. 

ARG 3. KLIST IS EITHER 1 OR 2 AND DECIDES WHICH SIDE OF 
KDATA YOU WANT THE DATA TO BE UNPACKED INTO 
ARG 4. KFLAG IS AN ERROR FLAG RAISED BY DATA. IF *0»N0 ERROR 
LOGICAL UNIT 1 IS USED BY DATA TO FIND BLOCK. 

EACH TIME DATA IS CALLED IT FINDS THE DESIRED BLOCK ON THE TAPE 
AND UNPACKS IT INTO EITHER KDAT A ( M » 1 ) OR KDATA(M,2> 

IN THE CALLING PROGRAM DIMENSION KDATA ( 4000 ♦ 2 ) » AND 
DECLARE KDATA COMMON. 

KDATA IS THE OUTPUT LIST AND IS REFERENCED BY KLIST 
SUBROUTINE DATA CALLS ON SR UNPACK AND SR FINDIT 
AFTER UNPACK THERE EXISTS ONE 1604 WORD/160 WORD 
ASSIGNMENT NOS USED IN DATA ,822 *823 , 824 ,825, 828 ,829 

824 FORMAT ( 48H1 ERROR OCCURRED IN FINDIT SUBROUTINE AT HEADER ) 

825 FORM AT (50X, 01 6) > 

828 FORMAT ( 33H0 I HAVE UNPACKED DATA HEADED BY ,016 ) 

830 FORMAT ( IX, 14, 3X,016) 

CALL FINDIT( IDENT, MAX, IFLAG) 

IF(IFLAG) 822,823,822 

822 PRINT 824 $ PRINT 825, IDENT $KFLAG=1$ GO TO 829 

823 JMAX=MAX/4+l 

CALL UNPACK ( IBLOCK,JMAX, KLIST ,KDATA( 1 , KLIST )) 

KFLAG = 0 $ PRINT 828 , IDENT 
DO 826 M= 1 , 4 

826 PRINT 830, M,KDATA(M, KLIST) 

J = MAX “ 3 

DO 827 M ■ J, MAX 

827 PRINT 830*M,KDATA(M, KLIST) 

829 CONTINUE 
END 
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SUBROUTINE FINOIT( IDENT*MAX* I FLAG) 

DIMENSION KDATA(4000»2) ♦ I BLOCK (1001 ) 

COMMON KDATA ♦ I BLOCK 

8 1 A FORMAT ( 53H1 A PARITY ERROR WAS DETECTED BUT RUN WAS NOT STOPPED) 

815 FORMAT (45H PARITY ERROR OCCURRED AT HEADER AS FOLLOWS ) 

816 FORMAT! 33X.016 ) 

818 FORMAT! 32H1 I HAVE LOCATED DATA HEADED BY *016 ) 

820 FORMAT !38H UNABLE TO LOCATE HEADING AS FOLLOWS ) 

MAXl c MAX/4 + 1 

ASSIGN 811 TO JUMP 

805 BUFFER IN ( 1 * 1 ) ( I BLOCK ! 1 ) * I BLOCK ( MAXI ) ) 

806 I F {UNIT ♦ 1 ) 806*807.808*810 

807 IF! IDENT-IBLOCK ( 1) ) 805 *8 13. 805 

808 GO TO JUMP* (811* 812) 

811 REWIND 1$ ASSIGN 812 TO JUMP $ GO TO 805 
810 IPAR=1 S GO TO 807 

812 REWIND 1$ GO TO 819 

809 PRINT 814$ PRINT 815$ PRINT 816*IDENT$ GO TO 817 

813 I F ! I PAR ) 809*817*809 

817 I FLAG * 0 $ PRINT 818*IDENT $ GO TO 821 

819 PRINT 820 
PRINT 816*1 DENT 
IFLAG-1 

ASSIGNMENT NOS USED HERE*805*807*808*809*810*8ll*812.813.814.815 

( 816*817*818.819*820*821 

821 CONTINUE 
END 
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TABLE E— li 
SUBROUTINE UNT *CK 





IDENT 




UNPACK 




/ / 
// 


ARCS I BLOCK ♦ JMAX ♦ KLIST»KDATA 




ENTRY 




UNPACK 




SUBS UNPACK CALLED BY SR DATA 


UNPACK 


SLJ 




*-* 






+ 




SIU 


1 


EXIT 










LIU 


1 


UNPACK 






+ JMAX=MAX/4+l » MOST* JMAX— 1 




LDA 


1 


0 






GETS ADDR OF CALLING ARGS 




SAL 




Ul 






+U1 IS JMAX ADDR 




ARS 




24 






I BLOCK PACKED IN DATA 1001 WDS 




INA 




-1 






+ 


U1 


SAU 




AADRS 






ADRS OF IBLOCK IN AADRS 




LDA 




** 






+FILLED IN BY 3 INSRUCTIONS BACK 




SAU 




BADRS 






STORE JMAX IN LOOP COUNT 




LDA 


1 


1 






+GO GET ADRS OF NEXT ARG 




SAU 




J4 






+STORE ADRS OF KDATA < 1 * KL I ST ) 




INA 




1 






ADVANCE ADRS 




SAL 




J3 






♦STORE ADRS OF KDAT A ( 2 *KLI ST ) 


• 


INA 




1 










SAU 




J2 






♦ DITTO FOR KDATAl 3 *KL 1ST ) 




INA 




1 










SAL 




Jl 






♦DITTO FOR KDATA < 4 *KL I ST ) 




INI 


1 


2 






SET CORRECT EXIT ADRS 




SIU 


1 


EXIT+1 






♦STORE IN EXIT INST 




SIL 


2 


EXIT 






SAVE INDEX 2 




ENI 


1 


2 


r 




♦FIRST DATA WORD IN IBLOCK 




ENI 


2 


0 




FIRST WORD IN KDATA IS DATA 


aadrs 


LDA 


1 


** 






+ IBLOCK ( J ) 




LRS 




12 






SHIFT 4TH WORD INTO A REG 


JI 


ORS 




36 






♦ RIGHT JUST *SIGN EXTEND 




STO 


2 


** 






STORE IN KDATA (I+4*KLIST) 




LRS 




12 






♦SHIFT 3RD WORD INTO Q 




ORS 




36 






RIGHT JUST ♦ S I GN EXTEND 


J2 


STQ 


2 


** 






♦STORE IN KDATA(I+3»KLIST) 




LRS 




12 






SHIFT 2ND WORD INTO A REG 


J3 


ORS 




36 






+ RIGHT JUST ♦ SIGN EXTEND 




STO 


2 


** 






STORE IN KDATA( 1 + 2* KL I ST) 


J4 


STA 


2 


** 






♦1ST WORD NOW RT.JUST +SIGN EXT 




INI 


2 


A 






1 = 1+4 


BADRS 


ISK 


1 


** 






+ISK ON JMAX 




SLJ 




AADRS 






J=J+1 REPEAT LOOP 


EXIT 


ENI 


1 


** 






♦RESTORE INDEX 1 




ENI 


2 


** 






RESTORE INDEX 2 




SLJ 




** 






+ JUMP OUT 
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IDENT 




UNPACK 


ARGS I BLOCK* JMAX »KL ISTtKDATA 




ENTRY 




UNPACK 


SUBR UNPACK CALLED BY SR DATA 


UNPACK 


SLJ 




** 


+ 




SIU 


1 


EXIT 






LIU 


1 


UNPACK 


+ JMAX “MAX/ 4+1 t MOST* JMAX-1 




LDA 


1 


0 


GETS ADDR OF CALLING ARGS 




SAL 




U1 


+U1 IS JMAX ADDR 




ARS 




24 


I BLOCK PACKED IN DATA 1001 WDS 




INA 




-1 


+ 


U1 


SAU 




AADRS 


ADRS OF IBLOCK IN AADRS 




LDA 




** 


+FILLED IN BY 3 INSRUCTIONS BACK 




SAU 




SADRS 


STORE JMAX IN LOOP COUNT 




LDA 


1 


1 


+GO GET ADRS OF NEXT ARG 




SAU 




J4 


+STORE ADRS OF KDATA < 1 ,KL I ST ) 




INA 




1 


ADVANCE ADRS 




SAL 




J3 


+STORE ADRS OF KDATA < 2 *KL I ST ) 




INA 




1 






SAU 




J2 


+DITTO FOR KDATA! 3*KLIST ) 




INA 




1 






SAL 




J1 


+DITTO FOR KDATA! 4 *KL 1ST ) 




INI 


1 


2 


SET CORRECT EXIT ADRS 




SIU 


1 


EXIT+1 


+STORE IN EXIT INST 




SIL 


2 


EXIT 


SAVE INDEX 2 




ENI 


1 


2 


+FIRST DATA WORD IN IBLOCK 




ENI 


2 


o 


FIRST WORD IN KDATA IS DATA 


AAORS 


LDA 


1 


** / 


+ I BLOCK ( J ) 




LRS 




12 


SHIFT 4TH WORD INTO A REG 


J1 


QRS 




36 


+RIGHT JUST *S IGN EXTEND 




STO 


2 


** 


STORE IN KDATA (I+4*KLIST> 




LRS 




12 


+SHIFT 3RD WORD INTO 0 




QRS 




36 


RIGHT JUST ,SIGN EXTEND 


J2 


STQ 


2 


** 


+STORE IN KDATA! I+3*KLIST) 




LRS 




12 


SHIFT 2ND WORD INTO A REG 


J3 


QRS 




36 


+ RIGHT JUST* SIGN EXTEND 




STQ 


2 


** 


STORE IN KDATA! I+2*KLIST) 


J4 


STA 


2 


** 


+1ST WORD NOW RT • JUST +SIGN EXT 




INI 


2 


4 


1 = 1+4 


SADRS 


ISK 


1 


** 


+ISK ON JMAX 




SLJ 




AADRS 


J=J+1 REPEAT LOOP 


EXIT 


ENI 


1 


** 


+RESTORE INDEX 1 




ENI 


2 


#* 


RESTORE INDEX 2 




SLJ 

END 

FINIS 




#* 


+ JUMP OUT 



-EXECUTE* 
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TABLE E-5 (Con't.) 



TUKEY SPEC 
XF ACT = .99626 A(0) 



TAU(SEC) 


AN( TAU) 


0 


1.00000 


.004 


. 46700 


.008 


-.36148 


.012 


-.48622 


.016 


. 00744 


.020 


.4101 1 


.024 


. 284c2 


.028 


-. 15434 


.032 


-.25875 


.036 


. 10725 


.040 


. 32453 


.044 


.07270 


.040 


-. 24842 


.052 


-.21618 


.056 


.13418 


.060 


.37584 


.064 


. 13704 


.068 


-. 26967 


.072 


-. 26584 


.076 


.08314 


.080 


.27917 


.084 


. 12765 


.088 


-. 16504 


.092 


-.21833 


.096 


.07653 


.100 


.28530 


.104 


.08008 


.108 


-. 19886 


. 112 


-. 15890 


.1 16 


.09799 


.120 


.23746 


.124 


.06567 


.128 


-.21386 


.132 


-. 17481 


.136 


. 12947 


. 1 40 


.22078 


.144 


.00913 


.148 


-. 17675 


.152 


-. 14032 


. 1 56 


. 06086 


.160 


. 17405 


. 1 64 


.01655 


. 1 68 


-. 16435 


.172 


-. 07287 


.176 


. 1 1040 


.180 


.12215 


. 1 84 


-.00197 


.188 


-. 12345 


.192 


-.09501 


. 196 


.09308 


.200 


. 16188 


.204 


-.00033 


.208 


-. 12095 


.212 


-.06221 


.216 


.04386 


.220 


.08959 


:Si8 


-:?W8 


.232 


-.07646 



V ESTIMATES 
.79463 



CPS ) 


XN(FREG) 


C 


.00749 


.ooc 


.01458 


.ooc 


' .01270 


.coc 


.01061 


.ooc 


.00873 


.coc 


.00686 


.ooc 


.00610 


.ooc 


.00614 


.ooc. 


.00617 


.ooc 


.00585 


.ooc 


.00691 


.ooc 


.00722 


.ooc 


.00484 


.ooc 


.00309 


.ooc 


.00344 


.ooc 


.00544 


.ooc 


.00642 


.ooc 


.00494 


.ooc 


.00499 


.coc 


.00670 


.ooc 


.00655 


.ooc 


.00568 


.ooc 


.00596 


.ooc 


.00661 


.ooc 


.00638 


.ooc 


.00557 


.ooc 


.00494 


.ooc 


.00591 


.coc 


.00934 


.000 


.01015 


.coc 


.01056 


.coc 


.01557 


.ooc 


.01536 


.ooc 


.01345 


.ooc 


.02101 


.coc 


.0256 1 


.000 


.02141 


.ooc 


.01787 


.coc 


.01 846 


.000 


.01716 


.00 0 


.01409 


.coc 


.01360 


.ooc 


.01429 


.000 


.01487 


.000 


.Cl 599 


.ooc 


.01616 


.coc 


.02042 


.ooc 


.0265 1 


.ooc 


.04243 


.ooc 


.07895 


.coc 


. 1 1 559 


.coc 


.09713 


.ooc 


.04315 


.ooc 


.01450 


.coc 


.00687 


.coc 


.00322 


.coc 


.00380 


.ooc 


.00677 


.ooc 


.00809 



T R U 

FREQ ( 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

i? 

58 
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.236 


.07577 


59 . 00 C 


.00541 


.240 


.07764 


60.00C 


.00286 


.244 


-.04312 


61 . GO C 


.00214 


.248 


-.09855 


62 .COC 


.00167 


.252 


-.0 390 8 


63.000 


.CO 147 


.256 


. 10244 


64 .COC 


.00157 


.260 


. 15790 


65 .COC 


.00228 


.264 


-. 000 1 2 


66 .GOC 


.00275. 


.268 


-. 1 4579 


67 .000 


. C02C5 


.272 


-.0607 1 


68 .00 0 


.00203 


.276 


.06367 


69. GOC 


.00241 


.260 


. 06875 


70.000 


. CO 1 90 


.284 


.01631 


71 .000 


.00147 


.288 


-.04 142 


72 .OOC 


.00139 


.292 


-. 03549 


73 .COC 


.00150 


.2 96 


.04668 


74 .OOC 


.CO 1 64 


.300 


.03850 


75.000 


.CO 1 54 


.304 


-.07314 


76. GOC 


.00127 


.308 


-.07468 


77. COC 


.00141 


.3 12 


.02519 


78 .OOG 


.00090 


.316 


.07875 


79. OOC 


.00323 


.320 


.07094 


80 .OOC 


.00973 


.324 


-.00860 


ei .COO 


.00816 


.328 


-. 10054 


82.000 


.C0151 


.332 


-.05049 


83.000 


.00019 


.336 


.04 195 


84 .COC 


.00035 


.340 


.01117 


85 .OOC 


.00019 


.344 


-. 03365 


86. OOC 


.00021 


.348 
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FIGURE E- 
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FIGURE E— 2 

SAMPLE GRAPHIC OUTPUT OF POWER SPECTRAL DENSITY PROGRAM 




LO 



\ 

ZP 

CO 



cc 



CO 

UJ 



CO 

tz 

UJ 

o 



cc 

cc 



o 

r iU 

2 2 Cl 
V C CO 

co co 
H I — 

Z.Z.OL 



CL 



CO CO 

c? 

c3 <m 

y y cc 

o <-> UJ 
‘P CO 
x v CO 



"O 



76 



APPENDIX F 



DEFINITION OF AUTO -CORRELATION FUNCTION 
AND POWER SPECTRAL DENSITY 



It is well known that under certain conditions a periodic 
function of a real variable can be represented by a Fourier series: 

f(<) = £ K C0S hUnA. V- b n SIN 

f) - ~ oO ' • 

where T is the period of the signal and the coefficients are given 
by: 

do = Tr * ; far r-Tf f m S'" 4i 

-T 

It follows from this that the power in a single periodic signal 

exists only at the frequency n/T cycles per second where n is any 

integer. The amplitude of the sinusoidal component of frequency 
2 2 

n/T is a^ + b^ . Thus the power in this component is 
2 2 

(a^ + b^ )/2. To express the power spectrum directly as a function 

of frequency, it is convenient to employ the Dirac delta function 
cf which has the following formal properties: 

d (iO -0 l) o ) ^ O CO% CO 0 



oO 

( J(lo- 0J o ) dw = f(<v e ) 



then the power spectrum can be expressed as: 






0 ) 



O- -<o 

The auto-correlation function of a periodic function is defined as: 






~r J fit) f(-6 + T)d.t 

-T 
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Computation of this quantity involves multiplication of the function 
by its time-shifted replica followed by integration over one period. 
Using the Fourier representation of equation (1) the correlation 
function of the periodic function f(t) is: 






R(rc )*T j 1 Z Un cos ~ + b n 

-Vl. na-ol — ' / _J 

^ Cos f brn S/N 

£(?) = "r t C" V* ^ COS ^ 



-Vi. 

£(T)=JT 3n* br£ c * s ZjrnV 

n- -oo ^ T 



U) 



The Fourier transform of a cosine is given by: 



JT(cos «J 0 i)= J COS co a t e~* ujt dt= ^(co-lo 0 ) (3) 

- Oq> 

thus, it is readily seen from equations 1, 2 and 3 that for a 
periodic function 



0(u>)* J R (T$ e'^'W'tr 

- CK> 

/in)-- Tw C?M e'^du 

-co 

that is, the power spectrum and auto-correlation function of a 
periodic signal constitute a transform pair. The equations are 
known as the Wiener -Khintchine relations. By letting't =0 we obtain 
directly the Parseval theorem 




This relation states that the power in a periodic signal can be 
obtained either by averaging the square of the signal over one period 
or by integrating the power spectrum [3]. 



78 



The definition of the auto -correlation function and power 
spectral density are quite clear for a periodic signal. However, 
a stationary random function is considered in this report. 

Let u(t) represent a random function of time. Normally a 
complete statistical description is necessary to define a random 
signal. For dynamic-load problems in aeroelast icity a time average 
concept may be used. The mean and the mean square of a function of 
time u(t) over a time interval 2T are: 




If the time average is independent of t^ and T when T is 
sufficiently large, then u(t) is said to be stationary. 

Assume that the function u(t) is observed for a long period 
2T, and that it vanishes outside the interval 2T. Further assume 
that 2T is greater than any characteristic time of the random 
function. Hence u^(t) represents all average statistical properties 
of the random function. 
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Since the truncation process assures absolute convergence and no 
finite discontinuities exist, the Fourier transform and inverse 
exist. Therefore, 



1 >) ! 



N 



(JO 



since u(t) is real 

where is the complex conjugate. 

The auto -correla tion function is defined as: 

„r 



( 5 ) 



TP fu(t) u (it't-)di (6) 

J-T 



H J c(oo Id* 



substituting the transform definition of u(t+t ) 

t - VP (vVr J *fr 

rearranging the order of terms,, 

*(r) • r? a '“‘ (Pm / Jdu> 

7 r *M = Ty(- W >- pr j <\M e '“' f «/t 



finally: 



define : 






It j J~(iv) Y^ , (w)e-‘'“ r d 



CO 



<pw- iz- -f\ T r H 

/?(?■) -i 



(?) 
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since by setting T =0 



#(o)= u l (i) = ± f <p(co)dco = f°0(io)cL^ (8) 

- oo A 



By definition g ^ - ->«.* £/* the imaginary 

part vanishes due to symmetry. 




By comparison with previous definitions 



<t>&) - -Jr j £(7>) tv T d U 



It is quite evident that the power spectral density is a 
Fourier cosine transform of the auto-correlation function and for 
a random function this is quite significant since the auto- 
correlation function can be easily measured and more clearly 
understood . 

Plotting (p( bj) against to produces a curve, the area under 
which represents the mean square of the random function as 
indicated in equation 8. The spectral character can appear either 
as a distributed form (type A) or as a dominant frequency (type B) 
or a combination of either. The figure below is in an approximate 
normalized form and illustrates the two types of responses. 
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